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The transverse mass spectra of O hyperons and <j) mesons measured recently by STAR Collabo¬ 
ration in Au-|-Au collisions at ,^snn = 130 GeV are described within a hydrodynamic model of the 
quark gluon plasma expansion and hadronization. The flow parameters at the plasma hadronization 
extracted by fitting these data are used to predict the transverse mass spectra of J/if and if' mesons. 


Recent measurements [1] of the energy dependence of 
pion and kaon production in central collisions of heavy 
nuclei (A-l-A) indicate [2] that the transient state of de- 
confined matter is created at the early stage of these col¬ 
lisions for energies higher than about 40 A-GeV, i.e. at 
high SPS and RHIC energies. Analysis of the hadron 
multiplicities measured in these collisions within statis¬ 
tical models at the SPS [3] and RHIC [4] energies shows 
that the chemical freeze-out takes place near the bound¬ 
ary between the quark-gluon and hadron phases. The 
values of the temperature parameter extracted from the 
data are similar for both energies, Th = 170 ± 10 MeV, 
and they are close to the value of the deconfinement tran¬ 
sition temperature at zero baryonic density estimated in 
the lattice QCD (see e.g. Ref. [5]). 

The analysis of the experimental data at the SPS [6,7] 
and a numerical modeling of the hadron cascade stage in 
A-l-A collisions at SPS and RHIC energies [8,9] indicate 
that the kinetic (i.e., particle spectra) freeze-out of the 
most abundant hadrons takes place at temperatures sig¬ 
nificantly lower than Th- Nevertheless, one expects that 
the kinetic freeze-out of some heavy and weakly inter¬ 
acting hadrons (e.g. Ul hyperons and (j), J/tf-i 4’' mesons) 
may occur directly at the quark-gluon plasma (QGP) 
hadronization stage or close to it. Thus for these hadrons 
the chemical and kinetic freeze-outs coincide and are de¬ 
termined by the features of the QGP hadronization. For 
Ul hyperons and (j) mesons this expectation is based on the 
results of “hydro QGP -b hadron cascade” approach [8,9]. 
For J/4’ and ip' mesons this is our suggestion [10-12], 
which is a straightforward consequence of the recently 
proposed statistical mechanism of charmonia production 
at the QGP hadronization [13-16]. Within this approach 
the transverse mass spectra of the H, J/ip and ip' mea¬ 
sured in Pb-t-Pb collisions at the SPS have been described 
successfully [11]. The transverse mass spectra of hyper¬ 
ons [17] and ip mesons [18] produced in central Au-bAu 
collisions at ^/snn = 130 GeV were recently measured 
by the STAR Gollaboration. These data allow to test 
our model in the new energy regime. They also give a 
unique opportunity to extract parameters of the QGP 
hadronization at RHIG energies and consequently pre¬ 
dict spectra of J/ip and ip' mesons. 

Within a hydrodynamical approach of the QGP 


hadronization the transverse mass spectrum of i-th 
hadron in the central rapidity region can be written as 
(see, e.g., Ref. [19]): 
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where y is the particle longitudinal rapidity and yxiO = 
tanh”^ vt is the fluid transverse rapidity. Rh and th 
are, respectively, the transverse system size and proper 
time at the hadronization (i.e., at the boundary between 
the mixed phase and hadron matter), ^ = r/Rn is 
a relative transverse coordinate. The particle degener¬ 
acy and fugacity are denoted as di and A^, respectively, 
rriT = \/Px T mf is the hadron transverse mass, Ki and 
Jo are the modified Bessel functions. Parameter in 
Eq. (1) (75 [20] for i = (p,Ul and 7^ [14,15] for i = J/ip, 
Ip') describes a possible deviation of strange and charm 
hadrons from complete chemical equilibrium (n^ = 2 for 
(p, J/ip, Ip' and n-i = 3 for fl). 

The spectrum (1) is obtained under the assumption 
that the hydrodynamic expansion is longitudinally boost 
invariant and that the freeze-out occurs at constant lon¬ 
gitudinal proper time r = {t is the time and 2: 

is the longitudinal coordinate), i.e. the freeze-out time 
t is independent of the transverse coordinate r. In order 
to complete Eq. (1) the functional form of the transverse 
rapidity distribution of hadronizing matter yxiO has to 
be given. A linear flow profile, yxiO = Vx^"^ ’ C; used 
in our model is justified by the numerical calculations of 
Ref. [9]. 

Thus, in our model, the QGP hadronization is de¬ 
scribed by the following parameters: temperature Th, 
“volume” thR/hj maximum flow rapidity fugacities 

\i, and saturation factors 7i. Note that the (p, J/ip, ip' 
have no conserved charges and Ai = 1 for these parti¬ 
cles. We use the fixed values of the parameters Th = 
170 MeV, 7s = 1.0, Xq~ = 1/Aq+ = 1.09 (note that 
Aq- = exp[(/rB — 3/is)/T], where ps and ps are, respec¬ 
tively, baryon and strange chemical potentials). These 
(average) values of the chemical freeze-out parameters 
have been found in the hadron gas analysis [4] of the full 
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set of the midrapidity particle number ratios measured in 
central Au+Au collisions at ^/snn = 130 GeV. The fit 
to the WT-spectra of hyperons [17] and (f) mesons [18] 
measured in central (14% for and 11% for </)) Au+Au 
collisions at ^/SNN = 130 GeV is shown in Fig. 1. The 
fit results are: = 0.74 ± 0.09, thRh = 275 ± 70 

fm^/c and y^/nd/ = 0.46. In the calculation of errors of 
the two free parameters of the model the uncertainties of 
Th (±5 MeV), 7 s (±0.05) and Aq- (±0.06) were taken 
into account. 
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FIG. 1. The hadron transverse mass spectra in Au+Au 
collisions at ^/snn = 130 GeV are shown. The points indi¬ 
cate experimental data for the Q [17] and 0 [18] measured by 
STAR. The model results are shown by full lines. 


A simple exponential approximation of the spectra is 
usually utilized to parameterize the experimental data: 

= C exp (- . (2) 

Note that in Refs. [10-12] an additional factor was 
present in the r.h.s. of Eq. (2). It led to smaller values of 
T* when fitting the same spectrum. The WT’-spectrum 
(1) may, however, deviate significantly from a purely ex¬ 
ponential one and its shape depends on the magnitude 
of the transverse flow and the mass of the particle. The 
normalization factors C and the inverse slope parameters 
T* in different intervals of rriT — m can be found from 
the (/), ft, J/ip and spectra given by Eq. (1) using the 
maximum likelihood method. The average values of T* 
for the rriT domains of “Iow-pt” (wt — m < 0.6 GeV) 
and “high-pr” (0.6 GeV< rriT — m <1.6 GeV), discussed 
in Refs. [9,12], are shown in Fig. 2. The values of T* ob¬ 
tained by fitting the J/'0 and ip' data in Pb+Pb 
collisions at 158 A-GeV (see Ref. [11]) are also shown for 
comparison. The observed increase of T* with increase 
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of the hadron mass is much stronger at RHIC than at 
SPS energies. It is caused by larger transverse flow ve¬ 
locity of hadronizing QGP at RHIG (wt — 0.44) than 
at SPS {vt — 0.19). The increase of T* is much more 
pronounced in “Iow-pt” region than in “high-pT” one. 
In our model the mr-spectra of charmonia are extraor¬ 
dinary affected by the stronger transverse flow at RHIC 
due to enormous masses of these hadrons. Thus, the data 
on J/ip and ip' production in Au+Au collisions, soon to 
be obtained at RHIC, should allow to test the hypothesis 
of their formation at the QGP hadronization. 



FIG. 2. The values of the inverse slope parameters T* for 
two different rriT domains - “low-py” (my — m < 0.6 GeV) 
and “high-pr” (0.6 GeV< rriT — m < 1.6 GeV) - in Au+Au 
collisions at ^^snn = 130 GeV are presented. They are found 
using Eq. (1) with Th = 170 MeV and ypp°‘^ = 0.74. For 
comparison, the values of T* extracted from fitting the data 
in Pb+Pb collisions at the SPS (Eq. (1) with Th = 170 MeV, 
yjnax _ Q 2 g^ gee Ref. [11]) are also shown. 

We note here that at present there exists an uncer¬ 
tainty in the estimates of the yc factor, therefore, the pre¬ 
dictions concerning charmonia multiplicities in Au+Au 
collisions at RHIC within statistical approaches signifi¬ 
cantly vary and their discussion goes beyond the scope 
of this letter. 

The “volume parameter” thR^ = ^{Th) extracted 
from the fit to the H and (p spectra defines the line 
Th = A{Th) ■ R/^ in the Rh-th plane. The allowed re¬ 
gion in the Rh~th plane can be estimated by varying the 
temperature parameter within its limits, Th = 165 MeV 
and Th = 175 MeV. The resulting lines are shown in 
Fig. 3. The transverse radius Rh = 5 + 7 fm and the 
proper time th =8 + 11 fm/c at the QGP hadronization 
can be estimated from the hydrodynamical calculations 
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of [9] for central Au+Au collisions at ^/s]^ = 130 GeV 
(see Fig. 3 in Ref. [9]). These model boundaries and their 
intersection with the Rh-th region found in our analysis 
are shown in Fig. 3. 
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FIG. 3. The lines th = A{Th) ■ of constant “volume 
parameter” A{Th) are shown: Th = 170 MeV corresponds 
to the dashed line, Th = 165 MeV and Th = 175 MeV cor¬ 
respond to the lower and upper solid lines, respectively. The 
dashed area is the intersection of the Rh-th region between 
the Th = 165 MeV and Th = 175 MeV lines with the region 
of Rh = 5 -r 7 fm and th = 8^11 fm/c estimated from 
Ref. [9]. 



^ loW—pT 

(MeV) 

rji* 

^ hiqh—pT 

(MeV) 

Refs. 

DATA p, p 

455 ± 105 

290 ± 40 

]23,25] 

Hydro-fRQMD 

480 

300 

[9J 

Single freeze-out 

315 

310 

[21,22] 

DATA A, A 

505 ± 60 

320 ± 30 

[24,26] 

Hydro-fRQMD 

440 

310 

[9] 

Single freeze-out 

360 

330 

[22] 


Table I. The values of inverse slope parameters T* for 
(anti)protons and (anti)lambdas in Au-|-Au collisions at 
y^sivjv = 130 GeV are presented. The experimental values 
are taken as the average ones over the STAR and PHENIX 
results (a difference in the results for particle and its anti¬ 
particle is small). 

Within our approach the mr-spectra of </>, fl, 
are assumed to be frozen at the space-time hyper-surface 
where the hadron phase starts. This assumption is justi¬ 
fied by the small hadronic cross sections and large masses 
of these particles (in addition, the WT-spectra of these 
hadrons are almost not affected by the resonance feed¬ 
ing). However, the mr-spectra of many other hadrons 
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FIG. 4. The points indicate the experimental mr-spectra 
of the A and A in central Au-I-Au collisions at 
y/SNN = 130 GeV measured by the STAR [24] and PHENIX 
[26] Collaborations. The “straight lines” are the exponen¬ 
tial approximations of the spectra with Eq. (2) in the low-pr 
(solid lines) and high-pT (dashed lines) regions. 


are expected to be significantly modified by hadronic 
rescattering. Contrary to this expectation it was recently 
postulated [21,22] that the simultaneous chemical and 
kinetic freeze-out in Au-|-Au collisions at RHIC occurs 
for all hadrons (a single freeze-out model). Do experi¬ 
mental data allow us to distinguish between these two 
approaches? 


The “hydro QGP -I- hadron cascade” approach [9] pre¬ 
dicts for central Au-fAu collisions at ^/SNN = 130 GeV 
that the hadron cascade stage modifies the mr-spectra 
of nucleons and A hyperons substantially. In particu¬ 
lar a large increase of the inverse slope parameter in 
the low-pT region is expected for these hadrons as a re¬ 
sult of hadronic rescattering and resonance decay effects. 
Thus measurements of (anti)proton and (anti)lambda 
TOT-spectra should allow to distinguish between the sin¬ 
gle freeze-out model and models which assume different 
kinetic freeze-out conditions for different hadrons. We 
performed the T* analysis of the present RHIG data from 
STAR [23,24] and PHENIX [25,26]. The resulting T* 
values are summarized in Table I together with the pre¬ 
dictions of the single freeze-out model [21,22] and “QGP 
hydro + hadron cascade” model [9]. The TOp-spectra of 
the A and A are also shown in Fig. 4. There are signif- 
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icant systematic differences between T* parameters ob¬ 
tained from the STAR and PHENIX data. In view of this 
fact, the values quoted in the Table I are calculated as an 
arithmetic average of both results, whereas the (system¬ 
atic) error was estimated to be a half of the difference 
between them. 

Despite the large uncertainties, the data seem to favor 
the “QGP hydro -I- hadron cascade” model over the single 
freeze-out model. Additional data in the low-p^ region 
and their theoretical analysis would be helpful to clarify 
presence of the hadron cascade stage and its influence on 
T*o^_pt of (anti)protons and (anti)lambdas. 

The results on rriT-spectra of charmonia in central 
Au-|-Au collisions at the RHIC energies are expected to 
be available soon. They should allow to test a statisti¬ 
cal approach to the charmonia production at the QGP 
hadronization in high energy nuclear collisions. In par¬ 
ticular, within this approach, we predict a strong (a 
few times) increase of the inverse slope parameter T* 
of the charmonia rriT-spectra at RHIG in comparison 
with that at SPS. The higher is the energy the larger 
inverse slope is expected due to increasing transverse 
flow of hadronizing QGP. Thus, at ^/snn = 200 GeV 
the increase of T* should become even more pronounced 
than at ^/snn = 130 GeV. Due to strong sensitivity of 
the charmonia spectra to the hadronization temperature 
and transverse flow velocity, their analysis should signif¬ 
icantly improve our estimate of these parameters. 
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